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ABSTRACT 

Using the WiFeS integral field spectrograph, we have undertaken a series of observations of 
star- forming galaxies in Compact Groups. In this first paper dedicated to the project, we present the 
analysis of the spiral galaxy NGC838, a member of the Hickson Compact Group 16, and of its galactic 
wind. Our observations reveal that the wind forms an asymmetric, bipolar, rotating structure, 
powered by a nuclear starburst. Emission line ratio diagnostics indicate that photoionization is the 
dominant excitation mechanism at the base of the wind. Mixing from slow shocks (up to 20%) 
increases further out along the outflow axis. The asymmetry of the wind is most likely caused by 
one of the two lobes of the wind bubble bursting out of its Hi envelope, as indicated by line ratios 
and radial velocity maps. The characteristics of this galactic wind suggest that it is caught early (a 
few Myr) in the wind evolution sequence. The wind is also quite different to the galactic wind in 
the partner galaxy NGC839 which contains a symmetric, shock-excited wind. Assuming that both 
galaxies have similar interaction histories, the two different winds must be a consequence of the 
intrinsic properties of NGC838 and NGC839 and their starbursts. 

Subject headings: galaxies: evolution, galaxies: individual (NGC838,NGC839), galaxies: interactions 



1. INTRODUCTION 

Galaxy evolution is influenced by a wide range of phe- 
nomena associated with galaxy interactions. These phe- 
nomena include the onset or quenching of star formation, 
the triggering of large scale galactic winds, gas strip- 
ping, tidal perturbations, or the activation and feed- 
ing of AGNs (Active Galactic Nuclei) and their asso- 
ciated feedback. Compact groups of galaxies (CGs) 
are dense environme nts prone to intense and repe ated 
galaxy interactions (jCoziol fc Plauchu-Fravnl 12007). In 
these groups, gas-rich, star forming galaxies evolve into 
gas-poor, quiescent objects. Galaxies may be first pro- 
cessed in CG-style environments, before being; trans- 
formed within filaments and larger clusters (Moran et al.l 
120071 : iCappellari et al.l 120111 : iMahajan et aljl2012h . The 
mid- infrared (mid-IR) colors of galaxies in CGs also sug- 
gest that galaxy evolution is ac celerated in CG eny iron- 
ments compared with the field (| Walker et al.ll2010[ ). 

The first catalog o f 100 low redshift CGs was compiled 
bv lHickso n (1982a,b) using the Palomar Observatory Sky 
Survey red prints. The gravitationally bounded nature of 
the majority of these Hickson Compact Groups (HCGs) 
was later confirmed using (1) spectroscopic observations 
(Hic kson et al.lfl992h and (2) X-ray observations detect- 
ing e mission from hot gas surroundin g the g r oup m em- 
bers (Ponma n et al.l IT996) . Recently. Ilovin o (2002) ex- 
tended the survey and compiled a list of 121 Southern 
Compact Groups (SCGs). CGs usually contain 5 to 15 
members, and have a radius of a few hundred kpc. 

CGs, and HCGs especially, have been subject to close 
scrutiny ever since their discovery. Radio ob servations 
of HCGs by iVerdes-Montenegro et al.l (|200lD revealed 
an overall deficiency of Hi gas of the order of 60%. 
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iMartinez-Badenes et al.l (j2Q12h confirmed this result and 
using CO lines detected an excess of molecular hy- 
drogen in spiral HC G galaxies of the order of 50%. 
IMartmez et al.l ((2010) found that HCGs contain a much 
higher fraction of low-luminosity AGN than the field or 
cluster environments. Yet, the overall star formation 
rate (SFR) of galaxies in HCG s is not significantly differ- 
ent f r om galaxies in the fie ld (Verdes-Montene gro et all 
1998: IMartmez et all 120101 : iBitsakis et al.l I2010L I2011D . 
Ide la Rosa et al.l (|2007l ) also observed that the HCG 
stellar population is on average older and more metal- 
poor than the field. These observations suggest that 
interactions-prone CGs quench star formation in the 
long term. The causal link between galaxy interac- 
tions a nd enhanced star formation activity is well estab- 
lished ([Barton et al.ll2Q00b [Bournaudl 1201 It I Alonso et al.l 
|2012[ ). but (1) these starburst phases are short lived, 
and (2) not all interac t ions t rigger a starburst episode 
(|Di Matteo et all I2007L 12008( 1. The competing role be- 
tween one-on-one interactions and the local environment 
of galaxies towards the type (and in tensity) of nuclear ac - 
tivity has been directly observed bv lSabater et al.l ([2013) 
in a statistical analys is of galaxies i n the Sloan Digital 
Sky Survey fSDSS: I^r¥^t^l2000h . 

In a so-called wet mergeiQ (e.g. Jos eph fc Wrightl 
Il985f h the galaxy's first encounter gives rise to tidal dis- 
turbances, often seen in the form of tidal tails and non- 
regular rotation curves. These tidal disturbances sub- 
sequently trigger a gas-infall towards the nuclei of the 
interacting galaxies. This strong and rapid infall in- 
duces shock waves that ionize the in-falling gas. The 
ionized gas cools, first via X-ray emission, later by op- 
tical and UV line emission, and finally by the excita- 
tion of molecular hydrogen formed in the warm, recom- 
bined post-shock gas. Molecular hydrogen excitation is 

2 for gas-rich galaxies, hydrodynamics is important, hence the 
term wet from hydro, the Greek root for water. 
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a more efficient cooling pathway than X-ray emission in 
certain systems. Stephanas Quintet is a prime example 
(jCluver et al.ll2010h . The extended formation of molecu- 
lar hydrogen (H2) creates large and dense clouds of warm 
H2. Shocks further induce gravitational instabilities in 
these gas clouds, which leads to strong episodes of mas- 
sive star formation. 

Theoretical models that incorporate the physics of the 
Interstellar Medium (ISM) over small scales (<1 pc) are 
lacking. Yet, they are critical to gain a detailed under- 
standing of the physics of galaxy evolution in e n viron - 
ments such as CGs. For example, iTeyssier et al.l ([20 10) 
present simulations with a resolution of 12 pc. They 
conclude that gas fragmentation into large clouds is the 
dominant star formation process. This mechanism can 
only be seen if the densest phases of the ISM are resolved. 

One of the main limitations of our current understand- 
ing of galaxy evolution is the type of instrumentation 
that has been available until recently. Galaxies are com- 
plex systems that need to be fully resolved both spa- 
tially and spectrally, and therefore observations with sin- 
gle fibers or with long-slit spectrographs are inadequate. 
Fortunately, Integral Field Spectrographs (IFS) now pro- 
vide astronomers with the optimum tool to tackle galaxy 
evolution at the sub-kpc level in the local Universe. The 
unintended discovery of a galactic wind during the first 
commissioning run of the Sydney- AAO Multi-object In- 
tegral Field spectrograph (SAMI: ICroom et al.ll2012f ) is 
one exa mple of the unique potential of IFS for galaxy 
surveys (|Fogartv et al.l 120121 ) . 

We are undertaking a series of IFS observations of star 
forming galaxies in CGs using the WiFeS instrument 
(jDopita et al.ll2007l ) on the Australian National Univer- 
sity (ANU) 2.3m telescope at Siding Spring Observatory. 
Because of its large field of view (25x38 square arcsec- 
onds), WiFeS is among the very few IFS currently ca- 
pable of observing the entire extent of CG galaxies in 
a reasonable amount of time. Our targets, which usu- 
ally span over 1 arcminute on sky, can be mapped com- 
pletely with just a few pointings. With this series of 
observations, we aim to investigate galaxy interactions, 
the formation of galactic winds and the onset of star 
formation episodes at kpc scales. This project is com- 
plementary to recent WiFeS observing ca mpaigns of iso- 
lated dwarf galaxies (Nic holls et al.ll2011 | . Nicholls et a l. 
2013, in preparation), strong mergers (Rich et al. 2 0lTh . 
LIRGs, (Luminous Infrared Galaxies; Ri ch et all 120121) 
and galaxies in cluster env ironments (jFarage et al.ll2010l 
l20l2UMerluzzi et al.ll2012l ). 

The low-redshift & high resolution approach of our 
project is complementary to other ongoing efforts to 
improve our understanding of galaxy evolution at high 
redshift, for example by using gravitational lensing to 
push the redshift limit of spatially resolved targets (e.g. 
iJones et al .1120101 1 Yuan et al.ll2011f ). Our study provides 
an important benchmark sample for these high redshift 
studies. We select ga laxies based on thei r mid-IR colors 
using WISE imag es dWright et all 12013) , UV color us- 
ing Galex images (M artin et al.l 120051 ). and optical color 
and morphology using Digital Sky Survey (DSS) images. 
This selection focuses the sample on galaxies with ongo- 
ing star formation. We use tidal features in the optical 
images to detect ongoing interactions. 

This series of observations is not designed nor intended 



as a survey. While we plan to target galaxies with a wide 
range of nearest-neighbour distances, located in a wide 
range of CGs sizes, our focus is on the detailed analysis 
of these galaxies and the physics at work over their whole 
spatial extent. 

In this article, we present the results of our IFS obser- 
vations of the galaxy NGC838 and reveal the true nature 
of its galactic wind. The paper is organized as follows. 
We introduce HCG16 and NGC838 in Section [2 Our 
observations and data reduction procedure are described 
in Section 03 Our method for extracting emission line 
fluxes and velocities is detailed in Section SJ Our re- 
sults are shown in Section [5] We discuss their impli- 
cations in Section [6] and summarise our conclusions in 
Section [71 In this paper, we assume Ho = 71 km s _1 
Mpc -1 , Vt M = 0.27 and fly =0.73, following the 7 years 
WMAP results (jLarson et al.l 1201 lh . 

2. NGC838 AND HCG16 

In the initial definition of Hickson (1982b), HCG16 was 
catalogued as containing 4 galaxies, with an on-sky diam- 
eter of 1.6 ar cminutes. Based on spec troscopic follow-up 
observations (jde Carvalho et aLlll997h . the total number 
of member galaxies has since increased to 7. Evidences 
suggest that HCG1 6 may be part of a larger structure 
(jRibeiro et al.l [l998h . where a core of 5 galaxies is sur- 
rounded by a halo of several others yet to fall deeper in 
the pote ntial well. HCG16 is one of the most active of the 
HCGs (jRibeiro et al.l[l996h . In Table HJ we list the dif- 
ferent group members and the ir activity type. Based on 
the H i structure of this group, IVerdes-Montenegro et al.l 
(2001) selected HCG16 as a typical example of the sec- 
ond phase of their 5-step evolution sequence for the Hi 
distribution in Compact Groups, where 30%-60% of the 
total Hi gas originally located in the halo of galaxies 
is dispersed in tidal features. In their radio map of 
the group, a clear Hi bridge is detected between the 
central 4 galaxies (NGC835,833,838,839) and NGC848. 
This bridge suggests that NGC848 may be an intruder 
galaxy that collided with HCG16 in a scenario similar 
to the archetypical Stefan's Quintet (a.k. a. HCG92, see 
iSulentic et aLlfeoQH lAppleton et al.ll2006D . 

NGC838 was part of the SINGG survey (Survey for 
Ionization in Neutral-Gas Galaxies. iMeurer et al.l [2006) 
consisting of Ha and R-band imaging for 468 galaxies se- 
lected in the HIPA SS survey (Hi Parkes All Sky Survey, 
iBarnes et al.ll200lh . The SINGG image of NGC838 and 
NGC839 is shown in Figure [TJ where the green dotted- 
line rectangle denotes the footprint of our WiFeS ob- 
servations. Both NGC838 and NGC839 have an associ- 
ated extended Ha emission region: a narrow and regu- 
lar structure for NGC839, and an asymmetric, extended 
structure for NGC838. The galactic wind nature of the 
Ha emission in NGC839, powered by an ongoing star- 
burst, was confirmed bv iRich et al.l (j2010f ). In Figure [2j 
we show the HST WFPC2 image of NGC838, and over- 
lay the footprint of our WiFeS observations as a dotted 
green rectangle. This image clearly shows a nuclear star- 
burst region surrounded by a dense system of filaments 
of dust and, presumably, molecular gas. At NGC838's 
distance of ~55 Mpc, 1" corresponds to 0.27 kpc. 

3. OBSERVATIONS AND DATA REDUCTION 
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TABLE 1 

Galaxy members of HCG16 and their activity 



Name 


Alt. name a 


R.A. b 


Dec. b 


cz b 


Activity 






(J2000) 


(J2000) 


(km s- 1 ) 




NGC835 


HCG16-a/l 


2 h 09 m 24.6 s 


-10°08'09" 


4073 


Seyfert 2, starburst 


NGC833 


HCG16-b/2 


2 h 09 m 20.8 s 


-10°07'59" 


3864 


LINER d , starburst 


NGC838 


HCG16-c/4 


2 h 09 m 38.5 s 


-10°08'48" 


3851 


starburst 


NGC839 


HCG16-d/5 


2 h 09 m 42.9 s 


-10°11'03" 


3874 


LINER d , starburst 


NGC848 


HCG16-3 


2 h 10 m 17.6 s 


-10°19'17" 


3989 


starburst 


KUG 0206-105 


HCG16-6 


2 h 09 m 06.0 s 


-10°19'13" 


3972 


starburst 




HCG16-10 


2 h 08 m 36.8 s 


-09°56'16" 


4027 


no emission lines detected 



a In lHicksonl Ill982b0 : l3e Carvalho etaT] (fl997h . The redshifts of HCG16-7,8,9,12,... rule out these galaxies as being members of HCG16. 
b R. A., Dec, and heliocen t ric velocity ex t racted from NED (Nasa Extragalactic Database). 
c See Rib eiro et al.l (fl996h : ITurner et all (I2001T) . 
d Low-Ionization Nuclear Emission Region 
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Fig. 1. — Ha (red in the online version) and R-band (blue 
in the online versio n) image of NG C838 and NGC839, from the 
SINGG survey (Meurer et al. 2006). The dotted square delineate 
our WiFeS observation footprint, covering the large majority of the 
Ha emission around NGC838. 



NGC838 was ob s erved with the WiFeS instrument 
(jDopita et all 120071 I2010D on the ANU 2.3m telescope 
at Siding Spring Observatory in October 2011. WiFeS 
is a dual-beam, image sheer integral field spectrograph 
with a spatial resolution of 1x0.5 square arcseconds and 
a spectral resolution of 3000 or 7000 depending on the 
chosen grating, with a field of view of 25x38 square arc- 
seconds. We have used the B3000 and R7000 gratings 
with the RT560 dichroic in our observations resulting in a 
complete wavelength coverage from ~3800A to ~7000A. 
Our blue spectra have a resolution of 3000 to cover all 
interesting emission lines blueward of 5 700 A and allow 
for line diagnostic analysis, while the red spectra have 
a resolution of 7000 to allow for a detailed study of the 
velocity structure of the Ka emission line. We used the 
1x2 binning mode, so that each spaxel in the final data 
cube covers a 1 x 1 square arcseconds area on-sky. 

Four WiFeS pointings were necessary to map NGC838 
and its extended Ha emission. A fifth field, centred on 
the galaxy, was observed to assist in the photometric 




Fig. 2. Wide-Field Planetary Camera 2 (WFPC2) image 
of NGC838 taken with the Hubble Space Telescope (HST, Rid.: 
10787, P.I.: J.Charlton). Numerous dense dust lanes can be easily 
identified. The dotted square represents the footprint of our WiFeS 
observations. This picture is a combination of 2 different images 
with different contrast settings showing both the bright core and 
the faint outskirts of the galaxy. The original images were obtained 
from the Hubble Legacy Archive. A color version of this figure is 
available in the online version of the article. 

calibration of the different fields used in the final mosaic. 
In Figure [3j we show the footprint of our WiFeS ob- 
servations superposed on a red band image of NGC838 
from the Second Digitized Sky Survey (DSS-2 p, along 
with t he NGC839 observations footprint of iRich et all 
([20T0h . Each field was observed 4 times, 1400 seconds 
each, resulting in 5600 seconds on target for each point- 
ing. A 700 second sky exposure was interleaved between 
each science exposure, while telluric and flux standard 
stars were observed throughout the night. The mosaic 
took 3.5 nights to complete, none of which was perfectly 
photometric. Seeing conditions ranged from 1.5" for the 
bottom left field, 1.8-2" for the central and bottom right 
field, and 2-2.5" for the top two fields. Data cubes ac- 
quired under better seeing conditions were not degraded 
to match the worst observed seeing of 2.5" . Each 4 expo- 

3 obtained from the European Southern Observatory Online 
Digitized Sky Survey Server 
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sures of a given field (shown in Figure |3j) were acquired 
close enough in time to ensure similar seeing conditions. 
Because we do a spaxel per spaxel analysis, subsequently 
degrading one field's spatial resolution to match the spa- 
tial resolution of the others is not required and would 
only result in data of poorer quality. 
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Fig. 3. Red band DSS-2 image of NGC838 (top right) and 
NGC839 (bottom left). The rectangles (in red in the online version) 
delineate our WiFeS observation footprint. The dashecMine rect- 
angles show the WiFeS observation footprints of Rich et al. (2010) 
for comparison. The DSS-2 image is copyright (c) 1993-2000 by 
the AAO and AURA. 

The data was reduced using the dedicated WiFeS 
IRAF pipel i ne (V ersion of 2011, December; see 
iDopita et al.l (j2Q10h for a complete description). Among 
some other minor updates, this version accounts for the 
fact that both the blue and the red WiFeS cameras are 
now read out through one amplifier only (instead of four, 
as previously). 

Each science exposure was reduced individually and 
converted into a bias subtracted, flat-fielded, flux and 
wavelength calibrated data cube. Bias subtraction 
was performed using a technique which accounts for 
the known temp oral variations of the WiFeS biases 
(|Rich et al.l 12010). Bias frames are acquired through- 
out the night. For a given science frame, a smooth plane 
is fitted to the closest acquired bias frame. Using these 
reconstructed bias frames instead of a standard master 
bias minimises the errors in bias subtraction. Telluric 
corrections were applied to the red spectra only. The 
blue and red spectra are joined at 5700A, with the blue 
spectra multiplied by a correction factor of the order of 
1.15 to account for a difference in the flux levels at this 
wavelength. 

We have used the IRAF task imcombine to assemble 
the final mosaic. The reduction pipeline produced 10 re- 
duced cubes (5 red and 5 blue) on the same wavelength 
grid. No further interpolation along the spectral axis 
was necessary when the 10 cubes were combined. We 
limited ourselves to integer spatial shifts, both for sim- 
plicity and due to the average seeing conditions during 
most of our observations. Because header coordinates 
cannot be trusted to the required accuracy of ~1", we 



implemented a custom Python program which uses the 
spatial flux correlations to best match the integrated Ha 
emission in the various overlap regions between fields. 
The final mosaic is 47 x 65 square arcseconds in size, 
and contains 3055 spectra. 

4. DATA ANALYSIS 
4.1. Stellar templates and emission line fitting 

The underlying gas physics in NGC838 can be under- 
stood using the emission line structures and intensities in 
the data cube. Extracting this information from each of 
the 3055 spectra contained within our final mosaic is con- 
ducted in a two step process. First, the underlying stel- 
lar continuum emission needs to be accurately removed. 
This step is especially critical if emission lines coincide 
with absorption features, which is the case for (but not 
only) the Ha and H/3 emission lines in NGC838. The 
second step requires a careful fitting of the sometimes 
complex line structures to accurately account for the to- 
tal line flux as well as to obtain an accurate description 
of the underlying gas kinematics. 

Performing this complex analysis on thousands of 
spectra is time consuming and requires an automated 
procedure. Here, we use the dedicated IDL pro- 
gram UHSPECFIT based on s tellar fitting routines by 
iMoustakas fc Kennicuttl (120061) . Th e cod e was origi- 
nally developed by lZahid et al.l (|201lD and iRupke et al.l 
(|2010f ). and later adapted by iRich et al.l (|2010l ) to be 
compatible with WiFeS data. UHSPECFIT is designed 
to handle the two different spectral resolutions in our 
data. The stellar contin uum is fitted using a combinatio n 
of stellar templates from lGonzalez Delgado et all ((2005). 
We are only interested in an accurate removal of the un- 
derlying stellar contribution, and therefore limited our- 
selves to th e solar metall i city, G eneva isochrones models, 
similarly to lRupke et al.l (|2010f ). A careful analysis of the 
underlying stellar population based on a careful stellar 
population model fitting would require higher S/N spec- 
tra, and is outside of the scope of this paper. 

The principal emission lines fitted are Ha, H/3, [O III 
AA4959,5007, [N n] A6583, [S n] AA6716,6731 and [O I 
A6300. Multiple gaussian components are fitted to the 
continuum-subtracted observed emission lines. The ve- 
locity and velocity dispersion of the different components 
are assumed to be identical for every emission line in a 
given spectrum. Our line fitting procedure can be sum- 
marized mathematically as follows : n emission lines for 
any given spaxel are described by n x m + m + m pa- 
rameters in total ; n x m intensities, m velocities and 
m velocity dispersions, where m G [1,2,3] is the number 
of component associated with a given spaxel. For each 
spaxel, gaussian components are classified based on their 
velocity dispersion, from the lowest (component 1, or cl) 
to the highest (component 3, or c3). 

If most spectra with S/N>5 display asymmetric line 
profiles, only a few have more than 2 components clearly 
identifiable. Each spectrum with S/N(Ha)>3 is fitted 
successively with a combination of 1, 2 and 3 gaussian 
components. We visually inspected every spectrum to 
determine the number of components that provides the 
best fit to each spectrum. Only 14 spaxels required 3 
gaussian components to provide an adequate fit to the 
emission line profiles. All other spectra could be well fit- 
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ted by one or a combination of two gaussian components. 

We implemented a simple f-test routine in IDL to en- 
sure that the decrease in the fit residual associated with 
an additional gaussian component is enough to .justify 
adding a component to the fit (e.g. IWestmoquette et al.l 
120071 ). We adopted a false-rejection probability of 0.05. 
In most cases, the f-test outcome was consistent with our 
visual selection, except for the spaxels located towards 
the center of NGC838, where the automated f-test indi- 
cated that three components may be justified. However, 
(1) two of these components have a width similar to our 
instrumental spectral resolution and are essentially unre- 
solved, and (2) fitting with only two components provides 
a fit residual within the noise level of the data in every 
case, so we favoured our visual selection over the f-test 
outcome in this region. In Figure HJ we show the number 
of fitted gaussian components used to describe the emis- 
sion lines structure for the 970 spaxels in our data where 
S/N(Ha)>3. The x and y-axis denote spaxels. A single 
spaxel corresponds to 1" or 0.27 kpc at the distance of 
NCG838. 
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Fig. 4. — Number of component required to best fit each spectra 
in the data cube with S/N(Hai)>3. Every spectrum was visually 
inspected to select the number of components best fitting its emis- 
sion line structures. 



In Figure [5j we show two fit examples that require 2 
(upper panel) and 3 (lower panel) gaussian components. 
In the Ha - [N II] panel, we also show the different gaus- 
sian components used to fit the Ha and N[ll] A6583 line 
profiles. 

The same number of components, different for every 
spaxel, has been used in the fits to all emission lines, blue 
or red. Because our spectra have both a high and a low 
resolution region (7000 in the red, compared to 3000 in 
the blue), multiple components visible in emission lines 
located in the high resolution regions are not always as 
clearly visible in the emission lines blueward of the stitch 
wavelength at 5700A. We have visually ensured that all 
red and blue emission line fit residuals are at the noise 
level of the data. However, individual gaussian compo- 
nents for the emission lines blueward of 5700A cannot 
be trusted because of their reduced S/N, and we will not 
use these individually. 



4.2. Galactic and extragalactic reddening 

Every spectrum in the final mosaic was corrected 
for g alactic extinction using the iSchlaflv fc Finkbeine"rl 
(|2011[ ) recalibration of the ISchlegel et all (j!998f ) ~ 
tinction map based on dust emission measured by 
COBE/ DIRBE and IRAS /ISSA. The recalibration as- 
sumes a IFitzpatrickl (jl999h reddening law with Ry =3.1 
and a different source spectrum than ISchlegel et al.l 
(fl998h . 

To estimate the extragalactic reddening along a given 
line of sight, we compute the total extinction in the V 
band Ay (in magnitudes) using the Ha to H/3 flux ratio. 
The resulting map is shown in Figure [6j where only the 
558 spaxels with S/N(Ha,H/3)>3 are displayed. 

NGC838 contains several large dust lanes clearly visi- 
ble in the HST image shown in Figure [21 In the reddening 
map, the impact of the most prominent dust lane, located 
to the South East of the galaxy center, can be readily 
identified as the region with Ay> 2.5 mag. The least 
obscured line-of sights are located South of this area, 
with Ay <1 mag. To the North West, the extinction is 
rather uniform with Ay =2 mag. This uniform extinc- 
tion region contrasts with the thick dust lanes visible to 
the North West of NGC838 in Figure [2 The lack of well 
defined dust lanes in this area of the reddening map in- 
dicates that the ionized gas is located above the central 
region of NGC838, consistent with an outflow scenario. 

The quoted Ay values in Figure [6] assume a case B 
recombination and associated intrinsinc Ha to H/3 flux 
ratio of 2.85, typical of Hll regions. This is a legitimate 
value because NGC838 is not hosting any AGN (which 
have been shown can in d uce larger intrinsic Ha to H/3 
ratios, see iKewlev et al.l (120061)). Ay va l ues h ave been 
computed using the Fischera fc Dopital (j2QQ5h extinc- 
tion curve with Ry=4.5, which provides self-consistent 
SFR estimates between U V, Ha and [Oil] luminosi- 
ties (jWijesinghe et al.ll2011l ). This theoretical extinction 
curve assumes the presence of a turbulent foreground 
dust screen, and (with Ry=4.5) is a good match to 
the emp i rical e xtinction curve for starburst galaxies of 
iCalzettil (j200l[ ). The exact equation of the e xtinction 
curve is deri ved i n the Appendix, see Eq. IA14I 

Using Eq. IA10| we corrected each spectrum for extra- 
galactic reddening using the associated Ha to H/3 flux 
ratio. The resulting emission line maps for the Ha, H/3, 
[O III] A5007, [N II] A6583 and [S II] A6716 + A6731 are 
shown in Figure [71 Only spaxels with S/N(Ha,H/3)>3 
for which an accurate reddening correction can be com- 
puted are shown. The corresponding flux density lower 
detection limit of emission lines in our observations is 
~lxl0 _1T erg s _1 cm -2 A -1 per spaxel. For a resolu- 
tion of R=7000, this line peak limit is equivalent to a 
minimal detectable flux (per spaxel) from an unresolved 
emission line (with a <25 km s _1 ) of ~1.3xl0 -17 erg 
s -1 cm~ 2 at 6000 A. 

The total Ha flux of NGC838 is \og(F Ha ) = 42.3 erg 
s _1 cm -2 after correcting for extragalactic reddening, 
and log^^Q,) = 41.5 erg s -1 cm" 2 before. This value 
is in perfect agreement with Me urer et al.l (|2006) who 
(correcting only for galactic extinction) obtain a total 
flux of \og(F Ha ) = 41.6 erg s" 1 cm" 2 for NGC838. 

5. RESULTS 
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Fig. 5. — Two examples of spectra (light-grey lines, in black in the online version) and their associated best-fit (black lines, in red in the 
online version). In the Ha - [N II] panel, the individual gaussian components used to fit the emission line are shown in dark-grey (blue in 
the online version). Respectively 2 and 3 components are required to fit the upper and lower spectra. Every column is labelled with the 
dominant emission line(s) it contains. The quoted wavelengths are the observed ones. 
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Fig. 6. — Reddening map of NGC838 in Ay magnitudes based 
on the Ha to H/3 line ratios for spaxels with S/N(Ho;,H/3)>3. The 
large dust lane seen in the HST image to the South East of NGC838 
is clearly visible. A color version of this Figure is available in the 
online version of the article. 

5.1. The underlying stellar population of NGC838 

In Figure [8j we show three spectra extracted from 
three apertures located in the center (top, thick line), 
to the West (middle line) and to the East (bottom line) 
of NGC838. The apertures are centred at the spaxel 
coordinates (23,34), (33,32) and (15,38) and are 4" in di- 
ameter. Each spectrum has been normalized at 4700A, 
and for clarity, the top and middle spectra have been off- 
set by +1 and +0.5 respectively. The wavelength scale is 
in the rest frame of the galaxy, assuming vngc838=3851 
km s _1 (see Tabled]). 

Despite a lower S/N in the lower two spectra, a clear 
signature of a young A-type stellar population is visi- 
ble in all three spectra. The strength of the absorption 
features associa ted with this y o ung ( ^500 Myr post- 
starburst, e.g. iBoehm-Vitensd 119921 : iStahler fc Pallal 



l2QQ5h stellar population is not diminishing away from 
the galaxy center. The spectrum extracted from the cen- 
tral region of NGC838 has a stronger blue continuum 
compared to the East and West ones. This signature of 
O/B-type star population is consistent with a very re- 
cent/ongoing starburst in the central region of NGC838. 
These young and hot stars are visible in the HST image 
of NGC838 as a central patch of blue emission. These 
features suggest that NGC838 has been subject to an 
episode of star formation throughout the whole extend 
of the galaxy some ~500 Myr ago, and that the activity 
has now settled down in the central regions of the galaxy. 

5.2. Velocity structure of the ionized gas 

The Ha flux, velocity and velocity dispersion maps for 
the cl (lowest velocity dispersion), c2 and c3 (highest ve- 
locity dispersion) groups of gaussian components fitted 
to the emission lines are shown in Figure [9l Velocities are 
in the NGC838 rest frame. A heliocentric velocity cor- 
rection of ~2 km s _1 was taken into account. Velocity 
dispersion have been corrected for instrumental resolu- 
tion (~lA at 6500A). 

Our extragalactic reddening correction relies on the H/3 
to Ha flux ratio, and only spaxels with S/N(H/3,Ha)>5 
can be corrected with enough accuracy to perform subse- 
quent line ratio analysis (see Sec. 15. 4ft . Consequently, our 
flux intensity maps contain less spaxels than the corre- 
sponding velocity and velocity dispersion maps, for which 
the cutoff criteria is S/N(Ha)>3. 

The cl and c2 velocity maps reveal a complex mixture 
of red- and blueshifted regions. A marked discontinu- 
ity in the cl flux map is visible at the transition region 
between 1 and 2 components. This discontinuity sug- 
gest that multiple components in the emission lines are 
present further away from the galaxy core than detected 
in our data, but cannot be resolved clearly due to the de- 
creasing S/N beyond ~10 spaxels from the center. The 
structure of the cl and c2 velocity and velocity disper- 
sion maps suggest a subdivision into four regions, delin- 
eated by dashed lines in Figure [9l We refer to them as 
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Fig. 7. Integrated flux maps for the H/3, [O III] A5007, Ha, [N II] A6583 and [S II] A6716 + A6731 lines after correcting for galactic and 
extragalactic extinction for spaxels with S/N(H/3,Hcn)>3. Our original flux detection limit is of the order of ~lxl0 — 17 erg s _1 cm -2 . A 
color version of this Figure is available in the online version of the article. 
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Fig. 8. — Spectra extracted from three apertures centered at the positions (23;34 - upper spectrum), (33,32 - middle spectrum) and (15,38 
- lower spectrum) and 4" in diameter. Each spectrum has been normalized at 4700 A. The upper and middle spectra have been offset by 
+1 and +0.5 respectively. Only 11 and 6 spectra with S/N(H/3,Hai)>3 where summed in the West and East apertures. A color version of 
this Figure is available in the online version. 



the South East Clump and the Northern, Central and 
Southern regions. 

The South East Clump — This region only comprises 
5 spaxels on the eastern edge of our mosaic. Its Ka 
line is well fitted by one gaussian component only. The 
associated velocity dispersion is low, of the order of 20-30 
km s _1 . A close inspection of the Ha SINGG image of 
NGC838 (see Figure [1]) reveals some possible features in 
this region, but the noise in the image makes these hard 
to identify with certainty. The clump is blueshifted, with 
a mean radial velocity equal to -122 ±6 km s _1 . 

The Central region — This area encompasses the center 
of NGC838. In the cl component, it has a low velocity 
dispersion (~25-50 km s _1 ) and a well defined signature 
of organised rotation. In the c2 component, a weaker 
rotation signature can also be identified, with a higher 
velocity dispersion of the order of 90 km s _1 . 

The Northern region — Most of this region is fitted with 
only one component, with the exception of a small region 
to the West. The associated velocity dispersion is of 
the order of 90 km s _1 , and a rotation signature can be 
identified in the cl velocity map. 

The Southern region — This region contains very high 
velocity dispersion values in both the cl and c2 maps, 
ranging from 90-150 km s _1 in the southernmost parts. 



This region does not display any sign of organised rota- 
tion in either the cl or c2 maps, but is strongly redshifted 
with velocities ~400 km s _1 , seen in the c2 velocity map. 
All spaxels fitted with three components are located in 
this region. 

In the next two Sections, we examine in detail the kine- 
matics of the Central, Northern and Southern regions. 

5.2.1. The rotation signature of the Northern and Central 

regions 

To best reveal the rotation signature of the Central 
and Northern region, we extract their associated rotation 
curve, and show them in the lower panel of Figure fTUl 
Each rotation curve is extracted along the relevant axis 
in the two top panels which show the Ha cl velocity map 
on the left and the stellar velocity map on the right. The 
stellar velocities are obtained from the stellar template 
fitting routine within UHSPECFIT. In each case, veloc- 
ity values are extracted in a 3 pixels band around each 
axis and averaged at every position. The associated 1-cr 
error is shown for every position. In the velocity maps 
(upper panels), the black points denote the location of 
the reference zero position on each axis, with the distance 
increasing positively towards the West (right). 

For the Central Region, we set the reference zero loca- 
tion at the spaxel coordinate (23;34). This corresponds 
to the Ha intensity peak location. The axis is inclined at 
95° West of North, our estimate of the P.A. of NGC838 



8 



Vogt et al. 



Flux [10~ 16 erg s~ 1 cm -2 ] 
50 100 150 200 250 



60 

50 

40 

30 

20 

10 


60 

50 

40 

30 

20 

10 


60 

50 

40 

30 

20 

10 



c1: 558 spaxels 
S/N(Hf3,Ha)>3 




c2: 437 spaxels 
S/N(H|3,Ha)>3 




c3: 15 spaxels 
S/N(H(3,Ha)>3 



Velocity [km s~ 1 ] 
-200-100 100 200 



d : 970 spaxels 
S/N(Ha)>3 




c2: 526 spaxels 
S/N(Ha)>3 



c3: 17 spaxels 
S/N(Ha)>3 



Vel. dispersion [km s 1 ] 
30 60 90 120 150 




_ c2: 526 spaxels 
S/N(Ha)>3 



_ c1 : 970 spaxels 
S/N(Ha)>3 



c3: 17 spaxels 
S/N(Ha)>3 



Northern 
region 

✓'Central region 



• Southern 
» region 

Sout East Clump 



10 20 30 40 



10 20 30 40 



10 20 30 40 



Fig. 9. — Ha flux (left), velocity (center) and velocity dispersion (right) maps of NGC838. Each row corresponds to a different gaussian 
component of the emission line, with maps labelled accordingly. Dashed lines mark the divisions between the four regions (Northern, 
Central, Southern and the South East Clump), labelled accordingly in the c3 velocity dispersion map (bottom right). The number of pixels 
in the maps and the associated S/N cut-offs are written in each panel. A color version of this Figure is available in the online version of 
the article. 



using our flux map s (see Figure [3. T his value is consis- 
tent with the NED (jHelou et al.lll991f ) values for the P.A. 
of NGC838 (81 degree in the r-band, 95 degrees in the 
Ks band). The (projected) rotation velocities range from 
+80 km s _1 to -140 km s _1 , spanning a total range of 220 
km s _1 . The central transition region between the red- 
and blueshifted side is very sharp, and the rotation curve 
flattens out rapidly further out. This well defined rota- 
tion signature most certainly arises from low-dispersion 
gas in orbit within the plane of the galaxy. To confirm 
this interpretation, we extract a velocity curve along the 
same axis using stellar velocities derived from the stellar 
velocity map (shown in Figure [TUJ upper right panel). 
This velocity map appears smoother than the Ha veloc- 
ity field in the central regions. The structure visible in 
the stellar velocity map to the East of the galaxy cen- 
ter is most likely a consequence of the small S/N in this 
area, which is also subject to high extinction values (see 
Figure [6]). The stellar rotation curve ranges from ^-100 



km s _1 to +80 km s _1 and is consistent with the gas ro- 
tation curve. Both rotation curves have an overall offset 
of -30 km s _1 with respect to the rest frame velocity of 
NGC838. 

The rotation curve in the Northern region is extracted 
along an axis inclined at 95° West from North and cen- 
tered at the spaxel coordinate (23,46). Its orientation 
is taken to be that of the axis along which the Central 
region rotation curve is extracted. The (projected) rota- 
tion velocity curve goes from ~-30 km s _1 to +110 km 
s _1 , spanning a total range of 140 km s _1 . The left- 
hand side of the curve is flat. The right-hand side, after 
peaking at 140 km s _1 , is decreasing towards a velocity 
of 100 km s _1 . The Northern region rotation curve has 
an overall redshift of +40 km s _1 with respect to the 
adopted rest frame velocity of NGC838, or +70 km s _1 
compared to the stellar rotation curve. Previous stud- 
ies of NGC838 which failed to distinguish between these 
multiple components led to erroneous estimates of its ra- 
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Fig. 10. — Top : cl velocity map of Ha (left) and stellar velocity- 
map (right). The A (full), B (dashed) and C (dotted) lines are the 
axes along which the radial velocity curves are extracted. Bottom 
: extracted radial velocity profiles. A color version of this Figure 
is available in the online version of the article. 



dial velocity. From our observations, it is clear that at 
least 2 distinct gas populations are present with very dis- 
tinct kinematic properties. From our measurement of the 
rotation curve of both the Ho; emitting gas and the stel- 
lar population, we find that the true radial velocity of 
NGC838 is most likely 30 km s _1 slower (-3820 km s" 1 ) 
than the current value found in the literature (—3850 
km s _1 in SDSS DR9, see lAhn et aHl2Q12h , highlighting 
the power of wide integral field spectroscopy for velocity 
field characterization. This radial velocity change does 
not challenge the membership of NGC838 to HCG16. We 
will discuss a possible explanation for the 70 km s _1 ve- 
locity offset between the Central and Northern regions 
in Sec. IO 

5.2.2. The redshifted Southern region 

To best reveal the complex velocity structure of the 
Southern region, we construct a position-velocity (PV) 
diagram of the cl, c2 and c3 velocity maps, shown in 
Figure [TTJ P V diagrams are a powerful tools which can 
simplify the visualisation of IFS 3-dimensional data by 
removing a spatial dimension. The concept is as follows. 
Every spaxel is projected onto a given axis, and assigned 
the corresponding distance from the reference location 
on this axis. This distance is then shown against the 
spaxel velocity in the PV diagram, which (with a care- 
fully chosen axis) can reveal complex structures in the 
data cube. 

The projection axis chosen to construct our PV dia- 
gram is shown in the two velocity maps in Figure [11] 
(upper panels). It is perpendicular to the radial veloc- 
ity axis in Figure [TUJ centered at the spaxel coordinate 
(23,34) and rotated 5° West of North. In the PV dia- 
gram (lower panel), red points correspond to the cl pix- 
els (lowest velocity dispersion), green crosses to the c2 
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Fig. 11. — PV diagram (lower panel) for the cl (lowest veloc- 
ity dispersion, black dots, in red in the online version), c2 (grey 
crosses, in green in the online version) and c3 (highest velocity dis- 
persion, black circles, in blue in the online version) velocity maps 
of NCG838. The projection axis is shown in the two velocity maps 
(upper panels), rotated 5° West from North and with the reference 
zero location at the position (23,34). 

The high redshift values present in the Southern region 
result in the large spread of c2 spaxels to the left of the 
PV diagram. The most redshifted spaxels have velocities 
slightly below 400 km s _1 . The distribution of both the 
cl and c2 spaxels does not follow any specific pattern for 
distances smaller than -5 pixels along the projection axis. 
This lack of structure is in strong contrast to the Cen- 
tral (-5 pixels to +7 pixels along the projection axis) and 
Northern (+7 pixels to +18 pixels along the projection 
axis) regions, which are organised in well defined struc- 
tures. The barred tubular structures are the (expected) 
result of the rotation of these two regions. The size of 
the envelope is equal to twice the (projected) rotation 
velocity of the gas. The offset between the Central and 
Northern region mean velocities of 70 km s _1 is reflected 
in the discontinuity visible at around +7 pixels along the 
projection axis in the red dot distribution. In addition, 
the smaller rotation velocity of the Northern region re- 
sults in a cylinder of red points thiner beyond a distance 
of +7 pixels from the reference point. In the Central re- 
gion, the green crosses trace another small barred tubu- 
lar structure, associated with the weak rotation of the c2 
component in this region seen previously in Figure [9l 

5.3. Velocity structure of the neutral gas 

Galactic winds can be detected and studied via ab- 
sorption lines, for example using the Na I AA5890,5896 
doublet (the Na D lines, see, e.g. Heckm an et al.l 2000: 
iRupke et al.l l2002f ) . Depending on the galaxy orienta- 
tion, Na D lines can probe either neutral gas within the 
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disk (in edge-on systems) or mat erial entrained w ithin 
the wind (in face-on systems, see Chen et al. 2010). The 
NED inclination for NGC838 based on Ks band obser- 
vations from the 2MASS survey is 46 degrees. This low 
inclination angle suggest that in NGC838, the Na D ab- 
sorption lines will primarily have their origin in material 
within the disk. 

We fit the Na D doublet for all spaxels with S/N(Na 
D)>5 (see Figure [5]). The resulting velocity map is shown 
in Figure [12] The Na D lines did not show any asymme- 
try or multiple line profile in any spaxel. They were 
always best fitted by single gaussian. For practical rea- 
sons, the Na D fits were performed separately from that 
of the emission lines using a dedicated Python routine. 
The quoted velocities are in the rest frame of NGC838, 
assuming vngc838=3851 km s _1 its radial velocity found 
in the literature. However, the colorbar is offset by 30 km 
s _1 (the offset we identified in Sec. 15.2. ip to allow for a 
better identification of inflows (red colours) and outflows 
(blue colours) with respect to the disk of NGC838. 
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Fig. 12. — Velocity map of the Na D doublet. The velocities are 
given in the rest frame of NGC838. The colorbar has been chosen 
so that it matches our measured radial velocity of NGC838. Black 
and white regions (red and blue in the online version) correspond 
to red- and blueshifted regions with respect to t he me an redshift 
of the stellar rotation curve of NGC838 (see Sec. 15. 2. If) . 



The velocity map reveals two distinct kinematic re- 
gions for the neutral gas with velocities of the order of 
+10 to -30 km s _1 . We do not detect any fast outflow sig- 
nature, consistent with the inclination of NGC838. To 
the East, a band of gas is redshifted by 20-30 km s _1 
with respect to the disk mean velocity. To the West, 
the gas velocity is of the order of the mean stellar veloc- 
ity of NGC838. A comparison with the HST image of 
NGC 838 and our extinction map in Figure [6] suggests 
that the East band may be the signature of gas inflow 
towards the galaxy center. This inflow is spatially co- 
herent with the large dust lane to the East of the galaxy 
center. This dust lane may be shielding the neutral gas 
from the main outflow. Further away from the dust lane 
(towards the North West), hints of a blueshifted neutral 
gas component suggest that the wind bubble's expansion 
is pushing the neutral gas away as it expands. This in- 
terpretation is subject to caution as the Na D velocity 



map is small. Higher S/N observations are required to 
confirm this scenario. 

5.4. Line ratio maps and diagnostics 

Specific line ratios can help understand the domi- 
nant excitation mecha ni sm of the outflow in NGC838 
dBaldwin etHI [19811; iVeilleux fc Osterbrockl 119871 : 
Veilleux & Rupke 2002; Matsubavashi et al. 2009). We 
use the [N n]/Ha, [S ll]/Ha, [O l]/Ha and [O lll]/Rf3 
ratios. These ratios are not strongly dependant on 
reddening corrections, a significant advantage in the 
case of NGC838 which is for the most part severely 
affected by intrinsic dust absorption (see Figure [6]). We 
show the line ratio maps for [N ll]/Ha and [S Il]/Ha in 
Figure [131 All three lines involved in these ratios are 
located in the red part of our spectra (high resolution), 
and the ratios for individual cl, c2 and c3 components 
can be computed. Because the c3 map is very scarce 
and does not contain critical information, we only 
show the line ratios for the cl (upper panels) and c2 
(lower panels) components. Only spaxels with S/N([N 
Il],Ha)>5 and S/NQS Il],Ha)>5 are shown. 
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Fig. 13. Line ratio maps for [N II] /Ha (left) and [S II] /Ha 
(right). The upper panels correspond to the cl component, and 
the lower panel to the c2 component. A color version of this Figure 
is available in the online version of this article. 

All four maps display the same general structure. Low 
line ratio values are located towards the center, and 
higher values towards the edges. For the cl component 
(lowest velocity dispersion), the map edges are thin (~2- 
5") with high line ratios. For the c2 maps, the high 
ratio value edges are thicker (~5-10"), but not located 
all around the maps. A tongue of low line ratios is ex- 
tending towards the South East all the way to the edge 



of the maps in both the [N ll]/Ha and [S II 
In Figure [I4j we show the [O l]/Ha and 



/Ha c2 maps. 
Om]/H/3 line 



ratio maps. Because (1) S/N([0 l])=5 and (2) the [O III] 
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and H/3 lines are located in the blue part of our spectra 
(low resolution), only the integrated, total flux ratios are 
computed. 
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Fig. 14. Line ratio maps for [O I] /Ha (left) and [O IIl]/H/3 
(right). A color version of this Figure is available in the online 
version of the article. 

The [O l]/Ha map shows a radial gradient with the 
lowest ratios centered on NGC838. We note the lack 
of any strong [O i] emission in the Southern region. 
For [O IIl]/H/3, the lowest ratios are located ~5" West 
of the Galaxy center. A core of small ratio values 
(log([0 III] /H/3) =-0.6) is surrounded by a thick envelope 
of higher values (log([0 III] /H/3) ^-0.3). 

We show the corresponding line ratio diagnostic dia- 
grams in Figure [15l Each data point corresponds to one 
spaxel in our data cube. Because any given spaxel has 
up to three different associated [N Il]/Ha and [S Il]/Ha, 
but only one value of [O III] /H/3, only the top two di- 
agrams have three data sets visible: cl (lowest veloc- 
ity dispersion, circles), c2 (diamonds) and c3 (highest 
velocity dispersion, triangle). The same [O III] /H/3 to- 
tal flux ratio is associated with all three (different) cl, 
c2 and c3 [N ll]/Ha and [S Il]/Ha ratios. In fact, no 
c3 point is shown because none of these spaxels have 
S/N([Q lll],H0)>5. T he different excitation regions de- 
fined in lKewlev et al.l (j2006f ) are shown. Each data point 
is coloured according to its velocity dispersion. The col- 
orbar extent and stretch is identical to Figure [9] For 
the lower diagnostic diagram, the velocity dispersion of 
each spaxel is taken as the Ha flux weighted mean ve- 
locity dispersion of the spaxels cl, c2 and c3 associated 
components. In the upper and central panel, a zoom 
on the region of interest is shown in a sub panel of 
the main diagrams. Photoionization mode ls computed 
with the St arburst99 dLeitherer et all 11999ft and MAP- 
PINGS III ([Sutherland fc Dopital 11993ft codes for ion- 
ization parameters Qg [6.50,6.75,7.00,7.25,7.50,7.75,8.00] 
are shown with stars symbols. The dotted line corre- 
sponds to solar metallicity, the full line to twice solar, 
and the dashed line to three times solar. 

All three diagnostic diagrams indicate that photoion- 
ization is the main excitation mechanism for the gas in 
the center of NGC838 (see Figure [H & ED • In the up- 
per and central panels, a clear dichotomy exists between 
low- and high- velocity dispersion spaxels. Two branches 
are visible in the upper panel - a low velocity dispersion 
branch (~30-50 km s _1 ) parallel to the Hll region delim- 
iting line, and a high velocity dispersion branch (~30-50 
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Fig. 15. — Standard line ratios diagnostic diagrams for NGC838. 
In the upper two panels, circles correspond to cl (low a) line ra- 
tios, diamonds to the c2 line ratios, and triangles to the c3 line 
ratios. In the lower panel, only the total flux line ratios are used, 
resulting in only one data point per spaxel. The Hl l, AGN and 
LINER region markers are from Kewlev et al. (2006). Points are 
colored according to their associated velocity dispersion. Stars are 
photoionization models for solar (dotted line), twice solar (full line) 
and three times solar (dashed line) metallicity and a varying ion- 
ization factor Q. Filled stars correspond to Q=8. A color version 
of this Figure is available in the online version of this article. 
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km s _1 ) extending towards and entering the composite 
region. We note that the former of these two branches 
contain in fact both low and high velocity dispersion 
points, while the second one only contains a high ve- 
locity dispersion component. These two branches can 
be interpreted as the respective signatures of (1) pure 
photoionization of gas within the disk of the galaxy and 
at the base of the wind and (2) mixing between pho- 
toionization and slow shocks along the wind axis. The 
photoionization branch is consistent with the twice so- 
lar metallicity photoionization model and an ionization 
parameter Q=7.0. In their study of ga laxy- wide shocks 
in late-stage mergers, Ric h et all ([20111 ) computed theo- 
retical line ratios for varying mixing values between pure 
photoionization and pure shock excitation. Comparing 
with Figure 10 in their paper, the composite line ratio 
branch that we detect in NGC838 is consistent with up 
to 20% shock fraction. 

In the lower panel, log([0 l]/Ha) is extending from - 
2.0 to -0.9, a wide range compared to the stretch of any 
other line ratio. Four data points also extend into the 
LINER region. These extreme ratio values are located 
away from the galaxy center (see Figure [T4|) . The in- 
creased [O i] line strength, reproducing to some extent 
a LINER signature, is a consequence of the geometry of 
the ionization bubble in NGC838. Towards the edge of 
the bubble, the geometrical depth along the ionization 
front is larger than towards the center. Overall, our line 
ratios are consistent with p revious measurements from 
Ide Carvalho fc Cozioll (fl999h for NGC838. 

6. DISCUSSION 
6.1. The global picture 

Previous, low-resolution observations of NGC838 sug- 
gested that it has a double nucleus with a 2" separa- 
tion and a complex velocity structure. These obser- 
vations wer e interpreted as the signature of an ongo- 
ing merger (jMendes de Oliveira et al1ll998[ ). We believe 
that Figure [2] and our WiFeS observations argue against 
this merger scenario, even if longer range interactions 
with neighbouring galaxies cannot be ruled out. We 
interpret the complex velocity structure detected pre- 
viously as the signature of a gaseous outflow powered 
by a nuclear star burst, in a sim ilar fashion to its neigh- 
bour NGC839 dRich et al.1120101) or the archetypical M82 
(jShopbell fc Bland-Hawthornl 119981 ). The nuclear star- 
burst is seen as an extended population of young, hot, 
blue stars in the HST image. We present in Figure [16] 
a toy model for the overall system which illustrates the 
different elements of our observations. 

From the lack of clearly defined features in the ex- 
tinction map (aside from the dominant dust lane to the 
South West, see Sec. 14. 2ft . we deduce that the gas re- 
sponsible for the Northern region emission seen in our 
observation is located above the galaxy plane, consis- 
tent with its outflow nature. Gas associated with this 
Northern region outflow has a well defined rotation (see 
Sec. I5.2.ip . The rotation is similar to that of the stellar 
component of NGC838, yet with a peak velocity 40 km 
s _1 slower and in a rest frame offset by 70 km s -1 . A 
galactic wind rotating slightly slower than its host galaxy 
is consistent with previous observations of other galac- 
tic winds (e.g. in M82, see lShopbell fc Bland-Hawthornl 
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Fig. 16. — Global schematic of the galactic wind in NGC838, 
with an expanding and rotating Northern wind bubble, a Southern 
wind bubble just entering its blowout phase, both of which are 
powered by a central starburst (not to scale). 



1998; Westmoquette et al. 2009), as well as with theoreti- 
cal models (jChattopadhv av" et al.ll2012[ ). As the wind ex- 
pands, conservation of angular momentum requires this 
to be the case. 

The HST image of NGC838 and its dust lanes makes 
the orientation of the galaxy un-ambiguous, placing the 
Southern side of the galaxy in front (see Figure [T6|). In 
that geometry, one would expect to see a blueshifted 
wind blowing towards the Earth in the Northern region. 
The 70 km s _1 mean redshift of the outflow in the North- 
ern region with respect to stellar disk mean velocity is 
therefore puzzling. We propose the following explana- 
tion for the wind mean redshift in this region. Thick 
dust lanes at the base of the wind act as barriers, and 
deflect the outflow toward the back of the galaxy. This 
scenario is motivated by the very thick dust lane to the 
South West of NGC838. It is also consistent with the 
velocity of the neutral gas in NGC838 (see Figure IT2]) . 
Neutral gas is falling- in towards the nucleus behind the 
dust lane, but transitions to a blueshifted velocity away 
from the dust lane (towards the N orth East), wh e re the 
wind bubble expands more easily. iCooper et al.l ((2008) 
have shown in their simulations that high gas (and dust) 
densities can significantly reduce the expansion of a wind 
bubble along the major axis of a galaxy. 

The Southern region contains highly redshifted, non- 
rotating gas, which is consistent with a leaking wind bub- 
ble or blow-out. In this scenario, the wind bubble, pow- 
ered by the nuclear starburst, has bursted out of the 
surroundin g Hi envelope into a lower density environ- 
ment (e.g. Veille ux et aTll2005D . The high pressure in- 
side the wind bubble is driving the rapid escape of the 
hot gas it contains and suppresses any existing wind ro- 
tation. We detected hints of multiple (> 3), distinct, 
velocity component in the Southern region, the possible 
signature of individual clumps being accelerated in the 
wind. Such filaments have been both detected in rea l 
winds, such as M82 (iShopbell fc Bland-Hawthornl 1 1998h 
and NGC839 ([Rich et al.1120101). or p redicted in theoret- 
ical simulations (|Cooper et al. 1120081 ). Their detection in 
the wind of NGC838 will require further observations 
with higher spatial resolution to be confirmed. We recall 
that the Southern region was observed under slightly bet- 
ter seeing conditions (1.5" -2"), which may explain why 
no hints of a filamentary structure is found in the North- 



The Galactic Winds of HCG16 



13 



ern region. 

In this toy model, the nature of the South East Clump 
described in Sec. 15.21 remains somewhat of a mystery. 
iCooper et al.l (j2008[ ) mention that in their simulations, 
the Ha filaments resulting from the breakup of gas clouds 
in the wind do not trace the true extent of the hot gas 
(defining the hot bubble). The South East Clump may 
therefore simply be the sign that the Southern wind bub- 
ble is more extended, and its breakout more generalised 
than shown in our toy model (see Figure [T6|). Another 
possible explanation, recalling that the signature of A- 
type stars is detected throughout the disk of NGC838, is 
that this Ha-bright clump is a remnant of another, ear- 
lier outflow episode related to a more global starburst. 
This over-density clump would become visible as it is 
photoionized by the central starburst in NGC838. Given 
its spatial location and blueshifted radial velocity of - 
122 ±6 km s _1 , the South East Clump does not seem 
to be related to the current Southern wind bubble and 
its highly redshifted material. A careful inspection of 
the SINGG Ha image of NGC838 reveals other knots of 
emission further to the East, outside of our WiFeS obser- 
vations. Follow-up observations are required to measure 
the radial velocity of these Ha knots, and, if confirmed, 
identify their nature (e.g. wind clumps or tidal debris 
fragments) and the origin of their eventual blueshift. 

What triggered the initial starburst in NGC838 re- 
mains to be clarified. At this point, the most likely 
explanation is that tidal perturbations resulting from 
the interactions of NCG838 with its neighbour NGC839 
triggered a gas infall fuelling the starburst. One could 
imagine anot her scenario (supported by the Hi map of 
HCG16, see iVerdes-Montenegro et all I2001D in which 
NGC848, colliding with HCG16 and passing close from 
NGC838 and NGC839, would trigger their respective nu- 
clear starburst and galactic wind formation. NGC848 is 
located ~14.2 arcminutes (204.5 kpc) to the South East 
of NGC838 (on-sky distance). Assuming an on-sky ve- 
locity of 400 km s _1 for the intruder galaxy (two times 
the typical HCGs velocity dispersion), the closest pas- 
sage would have occurred ~500 Gyr ago. This is in 
good agreement with the age of the A-type star popu- 
lation detected throughout the disk. The Northern wind 
bubble is extending roughly 5 kp c (on-sky) out of the 
galaxy plane, which according to Chattopadh vav et al.l 
(|2012l ) implies a wind age o f 37 Myr. Following the 
simulations of ICooper et al.l ([2008), the Southern wind 
bubble gas velocities lower than 400 km s _1 would indi- 
cate an even younger age for the wind, of the order of 
1-2 Myr. The ratio of the wind bubble radius to its ex- 
pansion velocity provides an alternative estimate of the 
wind expansion timescale. For the Northern bubble with 
a radius of ~5 kpc (on-sky) and an expansion velocity 
of ~100 km s _1 , we obtain an age of 49 Myr. For the 
Southern bubble with a radius of 7 kpc (on-sky) and an 
expansion velocity of 400 km s _1 , we obtain an age of 
17 Myr. These ages are only rough estimat es (see, e.g. 
IKoo fc McKeel Il992al lbl: IVeilleux et all 120051 and refer- 
ences therein for a more detailed treatment of the dynam- 
ics of galactic winds). There is a kn own delay (>10Myr, 
see lSharp fc B land-Hawt hornll2010l ) between a starburst 
episode and the onset of a galactic wind. Yet, the 2-3 or- 
ders of magnitude discrepancy between our first-order 



wind age estimates and the closest approach between 
NGC838 and NGC848 makes the latter a somewhat un- 
likely direct trigger for the activity seen in NGC838 to- 
day. It is plausible that a close encounter with NGC848 
triggered a global starburst episode in NGC838, which 
has since then settled down to the most central regions 
of the galaxy. The South East Clump may be one of the 
signature of this past interaction between NGC838 and 
NGC848. 

Line ratio diagnostics indicate that the base of the wind 
is photoionized by the central starburst, and that mixing 
with slow shocks increases further out along the outflow. 
Mixing is mostly detected in the Northern wind bubble, 
which is consistent with it bei ng still well contained with 
i ts sur rounding Hi envelope. Sharp fc Bland- Hawthorn! 
(|2010f ) noted in IFS observations of 10 low-redshift galac- 
tic winds that shock excitation usually dominates in in- 
stantaneous starburst driven wind. The largely pho- 
toionized starburst-driven wind of NGC838 therefore ar- 
gues in favour of our long-lasting, multi-phase starburst 
scenario (initially global, now central). 

6.2. NGC838 as a starburst-driven galactic wind 

Emission line splitting of a few hundred km s _1 
is a typical s ignature associated with many galactic 
winds (see, e.g. Heckman et al. 1990; Veilleux et al. 2005; 
Blan d- Hawthorn et al.l 120071 : [Sharp & Bland-Haw thornl 
20101 ). The classical interpretation is that the emitting 
gas is located on the surface of an open-ended bubble 
or conical structure. In NGC838, while we detect multi- 
ple kinematic components in the emission line profiles, a 
clear split is not evident in the Northern region. This re- 
sult indicates that the Northern wind bubble is still con- 
tained within its surroun envelope, and its interior not 
yet freely outflowing in the galaxy halo. NGC838 North- 
ern bubble is therefore caught early on in the galactic 
wind evolution sequence, a few Myr only after the win d 
starts to blow (jVeilleux et al.ll2005l : ICooper et alll2008f ). 

Large morphological differen ces exist b e tween 
different galactic winds (see IVeilleux et al.l [2005; 
IB land- Hawthorn et al.l I2007L for an extensive review) . 
These differences are related to the wind excitation 
mechanism (starburst or AGN), the host galaxy struc- 
ture, and the amount of gas present in the system. In 
the case of NGC838, the prime factor giving rise to the 
peculiar, asymmetric morphology of the outflow is that 
this wind is seen early on in its evolution. 

Limb-brightening is another feature seen in sev- 
eral galactic winds, for exam ple in NGC 307 9 
(jVeilleux et all fl99l. Arp 220 (iColina et al.l l200l . 
NGC 253 (IM atsubav ashi et al.ll2QQ9h or the Milky Way 
(jBland-Hawthorn fc Cohenl I2003D . In NGC838, we do 
not see any clear evidence for limb-brightening, sug- 
gesting that the wind bubbles are not hollow, but 
rather filled with the emitting gas. It is understood 
that optical emission in galactic winds originate from 
disk material entrained in the outflow, an interpreta- 
tion c onfirm ed by theoretical simulations (Co oper et al.l 
2008, l2QQ9h . The clear rotation signature that we de- 
tect in the Northern wind bubble of NGC838 is con- 
sistent with this entrained origin of the emitting gas. 
We note here that the wind in NGC838 is most cer- 
tainly very dissimila r to the limb-darkened wind in 
the Circinius galaxy (Vei lleux fc B land-H awthornl 1997: 
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ISharp fc Bland-Hawthornl I2010D . which is a Seyfert 2 
galaxy with a much more evolved wind. 

6.3. Comparison with X-ray observations 

X-ray emission from HCG16 has been reported by 
many authors in the past (e.g. iSaracco & Ciliesi 1995; 
Pon man et all 119961 : Turner et al.l [20 01: Bel sole et all 
l2003f ). Most recently, besjardins et al l (|2012h presented 
a new analysis of arc hival Chandra X-ray Observatory 
(jWeisskoof et al.ll2002h observations of HCG16 (ObsID: 
923, P.I.: Mamon) as part of a larger study of diffuse 
X-ray emission in CGs. They provide emission maps of 
the diffuse soft X-ray emission (0.2-2.0 keV) in HCG16 
which is clearly concentrated around the galaxy members 
of the group. The spatial correlation between the soft X- 
ray emission and the Ha emission is evident for both 
NGC838 and NGC839 in the color mosaic of the group 
(combining R-ba nd, Ha and soft X-ray images). Keep- 
ing in mind that [Pes jar di mTet al.l (|2Q12h have smoothed 
their Chandra data, we note that the soft X-ray emission 
associated with NGC838 is surrounding the Ha emis- 
sion by ^5-10" in all directions except to the South-East, 
where the X-ray emission vanishes more rapidly. 

This lack of soft X-ray emission is consistent with the 
leaking Southern wind bubble scenario described pre- 
viously. The low density and high temperature of the 
gas comp osing galactic winds is hard t o detect directly 
(see, e.g. IB land- Hawthorn et al.l 120071 and references 
therein). Instead, galactic winds are usually seen indi- 
rectly through photoionized, shock-excited or absorbing 
neutral material entrained from the disc. In NGC838, 
the wind in the Northern bubble has not yet burst out of 
the surrounding gaseous envelope. The wind's encounter 
with the surrounding material at the bubble's edge drives 
1) a forward shock wave in the surrounding gaseous ma- 
terial, and 2) a reverse shock wave in the galactic wind. 
With velocities of the order of several hundred km s _1 
(i.e. of the order of the outflow velocity at the base of 
the wind), the reverse shock wave excites the wind ma- 
terial to X-ray temperatures which then cools via soft 
X-ray emission. The shocks will be brightest when the 
galactic wind is pushing through a denser surrounding 
environment, as the plasma emissivity is a function of 
the plasma density (the energy lost via two-body interac- 
tions _scjy£s_aj^ density squared. See, 
e.g. JDopita fc Sutherland! 120031 . p.143-144). If the wind 
material can leak into a low-density region of the sur- 
rounding medium, as is likely the case for the Southern 
wind bubble of NGC838, then the shock- heated X-ray 
plasma would rapidly expand. This expansion will lead 
to a decrease in the density and emissivity of the plasma, 
that would then become too faint to be detected. 

6.4. NGC838 versus NGC839 

The circumstances required to trigger galactic winds 
are not yet fully understood. Identifying the necessary 
ingredients and actions resulting in a large scale galactic 
wind is complicated by the fact that both the galaxies 
content and their environment are likely to play a role 
in the process. From that point of view, the galactic 
winds in NGC838 and NGC839 may represent a unique 
opportunity to improve our understanding of the forma- 
tion mechanisms of galactic winds as a function of the 



galaxy characteristics, su ch as the disk thick ness or the 
location of the starburst (jCooper et al.ll2008l ). 

Located in the same region of the same CG, NGC838 
and NGC839 are subject to very similar environments 
and interaction histories (including with the intruder 
galaxy NGC848). The on-sky disposition of galaxies 
within HC G16 and its Hi distribution ( see the VLA 
Hi map of IVerdes-Montenegro et aD 12001) suggest that 
NGC838 and NGC839 are likely to be their respec- 
tive closest interacting partner. Any difference in the 
galactic wind located in these galaxies must therefore 
be linked to intrinsic differences in the galaxies them- 
selves. The asymmetric, photoionized wind of NGC838 
is in strong contra st with the sym metric, shock-excited 
wind in NGC839 (jRich et al.l 120101 . see also Figure [Q). 
Furthermore, the conical morpholog y with constant out- 
flow velocity inferred by iRich et all ([2010) for the wind 
in NGC839 is very different from the wind morphol- 
ogy in NGC838. Understanding the differences between 
NGC838 and NGC839 is therefore key to understanding 
the differences in their galactic winds. 

For both galaxies, resolving the star formation history 
and identifying the underlying stellar population needs 
to be addressed carefully. For example, it may be that 
NGC838 and NGC839 have experienced episodic, out- 
of-phase starburst events giving rise to the differences in 
their galactic winds. Confirming or disproving this sce- 
nario requires high S/N measurements of the stellar con- 
tinuum. High resolution Hi maps would also be beneficial 
to understand where the gas reserve is located and what 
is the exact extent of the interaction between NGC838 
and NGC839. 

7. CONCLUSION 

We presented our IFS observation of NGC838 and its 
galactic wind. These are the first results of a series 
of observations targeting star forming galaxies in Com- 
pact Groups with the WiFeS instrument on the ANU 
2.3m telescope at Siding Spring Observatory. Our ob- 
servations reveal the complex signature of an asymmet- 
ric galactic wind in NGC838. Emission line ratio maps 
and diagnostic diagrams show that photoionization is the 
main excitation mechanism at the base of the wind, with 
mixing from slow shocks (up to 20%) increasing away 
from the galaxy center along the outflow axis. The free 
flowing gas in an open-ended Southern bubble compared 
to the contained and rotating gas in a closed Northern 
bubble gives rise to the asymmetry of the wind. The 
closed Northern bubble is a strong indicator that the 
Northern wind is caught early (a few Myr) in the galac- 
tic wind evolution sequence. We also find kinematic evi- 
dence supporting a scenario where the Northern outflow 
is blocked and/or redirected at its base by thick dust 
lanes present in NGC838. 

The presence of A-type stars throughout NGC838 sug- 
gests that the galaxy has been subject to a global episode 
of star formation some 500 Myr ago. The photoion- 
ized nature of the wind is consistent with this non- 
instantaneous starburst scenario, which started glob- 
ally and has now settled down in the center of the 
galaxy. The Hi map of HCG16 suggests that the intruder 
galaxy NGC848 may have triggered the original galaxy- 
wide starburst, as it collided with HCG16 and flew by 
NGC838. 
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Finally, the photoionized, expanding, rotating, asym- 
metric wind of NGC838 is in strong contrast with the 
symmetric and shock-excited wind of the neighbouring 
galaxy NGC839. Because they are subject to the same 
environment, and (most likely) the same interaction his- 
tory, the differences in these two galactic winds are 
mostly due to intrinsic differences in their host galaxies. 
NGC838 and NGC839 therefore represent a unique pair 
of galaxies holding critical clues regarding the formation 
mechanisms of galactic winds. 
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APPENDIX 
DUST REDDENING CORRECTIONS 

There exist many ways to correct for extragalactic reddening, which is in itself not very well defined. Unfortunately, 
the literature is rather confusing, and it is often left to the reader to reconstruct by himself the exact solutions 
implemented in a given paper. Many authors also confuse extinction, which applies along a single sight line, with 
mean attenuation, which is a s tatistical sum of the attenuatio ns produced by a foreground dust screen with a fractal 
distribution (see the papers bv lFischera fc Dopitall2005l I2Q11L for the mathematic al treatment of this ). 

Here, we bas e our extragalactic reddening corrections on the work of Fisch era fc Dopital ([2005), which 
Wijesin ghe et al.l (1201 If) have sho wn provides the best agreement between star formation rate indicators in 
GAMA galaxies (Driver et al. 20 091). The mode l assumes that the attenuation is caused by a distant isothermal and 
turbulent dust screen. [Fischera fc Dopital (j2005[ ) provide several theoretical estimates of the relative attenuation as a 
function of wavelength which allow to correct a given spectrum based on its associated Ra/K/3 ratio. 

For clarity, let us deriv e the correction eq uation applied to our observations, starting from first principles. Let 
us adopt the formalism of lOsterbrockl ([1989) (see their Chater 7). Given a luminous source subject to attenuation, 
we define F\ the observed flux at a given wavelength A (in the source reference frame), and F\,o the theoretical, 
un-reddened flux. Then, 

^A = ^A,o-e" r \ (Al) 

where t\ is the opacity along the object line-of-sight. Both scattering and abso rption of the object light by the 
intervening dust are likely to play a role, and in fact, t\ = t\ a 5 S + t\ sca (see, e.g. iDopita fc Suther land 2003, p.296). 
With m\ and the observed and theoretical, un-reddened magnitude of the object, the extinction on the line-of- 
sight A\ for a given wavelength (in magnitude) is 

A x = m x - m A ,o = 2.5 log(e)r A . (A2) 

The actual shape of the opacity curve is highly debated. It is generally accepted as a simplification and based on 
galactic observations that t\ can be written as 

rx = cfx, (A3) 

where the shape of the curve is described by a universal function fx, and c is a multiplicative constant varying for 
every line-of-sight. Determining c can be done by comparing the observe d rat io of two given lines to the theoretically 
expected value, for example the Ra and H/3 lines. Using Eq. lAll and Eq. IA3| we have 



FHa = F Hai0 # e - c (f Ha -f Hfi ) 
Fh/3 Fh/3,0 



(A4) 
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and re- arranging the terms, 

c = -(f Ha ~ fnp)- 1 ■ (log(e))- 1 • log ( f Ha/ / l HP ) ■ (A5) 

Assuming Case B recombination, Fhcx,q/Fh/3,q = 2.85. However, it has also been sho wn that regions close from active 
galactic nuclei may have a higher intrinsic ratio value of 3.1 (e.g. iKewlev et aI1l2006[ ). 

To find de-reddened flux value, we substitute the above expression for c in Eq. IA1| which with simple algebra, and 
after re-arranging the terms slightly, becomes 



(A6) 



k Fhoc,o/Fh(3,0, 

where fx is the only unknown. Following Fisc hera fc Dopital (|2005f h we can write 

E x -v = A x - Ay 
Eb-V Ab — Ay ' 

where Ex-y is usually refereed to as the color excess between two bands, and Ex-y / Eb-v is the relative color excess. 
Introducing Ry = Ay / 'Eb-v, Eq. I A7I becomes 

Ex-v 



(A7) 



E B - 



-V 



Ry ■ Eb-v- (A8) 



Therefore, 



Finally, substituting Eq. IA9I in Eq. IA61 we find the reddening correction function, 



JX _ T X _ A X _ E B - V + U V 



fHa ~ fHf3 T Ha T H p A Ha - A H p E^LzZ _ Em^v 

h/B — V h/B-V 



FHa I Fh/3 



-v 



-E7 



F x ,o =F x -( - ""'J:* Eb ~ v " E *~ v • (A10) 

\-FHa, 0/^^/3,0 / 

A value of Ry = 4.3 results in an attenuation curve very simil ar to that defined empirically by iCalzettil (|2001f ) for 
starburst galaxies (jFischera et al. 2003; Fisc hera fc Do pita 2005). We therefore adopt the relative extinction curve 
with Ry =4.5 and Ay = 1 of Fisch era fc Dopital (j2005l ). From their Table 1, it is possible to extract the value of 
Ex-v I Eb-v for various values of A. From the few data points they provide in their paper, we derived a 4 th order 
polynomial (using a least-square fitting routine) allowing us to obtain an accurate estimate of Ex-v / 'Eb-v f° r an y 
value of A g[2480A ; 12390 A1. The derived function is shown in Figure [T71 (black line) alongside with the points from 
iFischera fc Dopital ([2005), and can be written as 

f A ~ y ^ -4.61777 + 1.41612 • A" 1 + 1.52077 • A" 2 - 0.63269 • A" 3 + 0.07386 • A" 4 (All) 

Eb-v 

with A in /am. We can then write Eq. IA10I as 

} 0.76- ( 5^7+4.5) 

K FHa,o/FH/3,0 

As a final remark, we note that w e ca n also write the expression allowing to convert an flux ratio to an Ay 

magnitude, by substituting c in Eq. IA2| which gives 

A V = -2.5 log ( J; Ha % HP ) • 

\^Ha, Q/ £H(3,Q J J Ha ~ jHfi 

Av = -2.5 log ( Jj^IlEL_) . ^ 8.55 log ( I Ha, ^ p ) (A14) 

tLB-V tuB — V 
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